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ABSTRACT 

An  experimental  investigation  of  the  response  of  a hypersonic 
turbulent  boundary  layer  to  a step  change  in  surface  roughness  has 
been  performed.  The  boundary  layer  on  a flat  nozzle  wall  of  a Mach 
6 wind  tunnel  was  subjected  to  abrupt  changes  in  surface  roughness 
and  its  adjustment  to  the  new  surface  conditions  was  examined.  Both 
mean  and  fluctuating  flow  properties  were  acquired  for  smooth-to- 
rough  and  rough- to- smooth  surface  configurations. 

The  boundary  layer  was  found  to  respond  gradually  and  to 
attain  new  equilibrium  profiles,  for  both  the  mean  and  the  fluctuating 
properties,  some  10  to  256  downstream  of  the  step  change.  Mean 
flow  self- similarity  was  the  first  to  establish  itself,  followed  by  the 
mass  flux  fluctuations,  followed  in  turn  by  the  total  temperature 
fluctuations. 

Use  of  a modified  Van  Driest  transformation  resulted  in  good 
correlation  of  smooth  and  rough  wall  data  in  the  form  of  the  incom- 
pressible law  of  the  wall.  This  is  true  even  in  the  nonequilibrium 
vicinity  of  the  step  for  small  roughness  heights. 

The  present  data  are  found  to  correlate  well  with  previously 
published  roughness  effect  data  from  low  and  high  speed  flows  when 
the  roughnesses  are  characterized  by  an  equivalent  sand  grain  rough- 
ness height. 

Existing  correlations  based  on  low  speed  data  were  found  to 
be  unsuccessful  in  predicting  the  effect  of  this  roughness  on  the  skin 
friction  and  velocity  profile.  The  indiscriminate  use  of  low  speed 


iii 

roughness  effects  correlations  to  predict  the  effects  of  roughness 
on  supersonic  and  hypersonic  flows  must  therefore  be  regarded  as 
a procedure  subject  to  gross  errors. 

Significant  pressure  and  temperature  history  effects  were 
observed  throughout  the  boundary  layer.  The  existence  of  these 
effects  was  found  to  create  a nozzle  wall  boundary  layer  whose 
properties  were  far  different  than  those  in  a boundary  layer  on  a flat 
plate  in  the  freestream,  raising  questions  about  the  validity  of  simu- 
lating the  flat  plate  boundary  layer  with  the  nozzle  wall  boundary 
layer. 
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INTRODUCTION 

The  effect  of  surface  roughness  on  the  cliaracteristics  of  flow 
over  a surface  has  long  been  of  interest,  especially  to  tiiose  people 
who  are  involved  in  the  design  of  veliicles  which  operate  in  or  on  the 
water  or  within  the  ^ttrnosphere.  It  is  well  known  that  the  presence  of 
surface  roughness  can  significantly  alter  the  drag  and  heat  transfer 
characteristics  of  a surface  and  can  even  cause  considerable  modi- 
fication of  the  flow-field  about  a body,  compared  to  the  smooth  wall 
case.  Quantitative  knowledge  of  these  roughness  effects,  and  insight 
into  the  physical  phenomena  which  give  rise  to  these  effects  is  essen- 
tial for  the  realistic  design  of  new  vehicles,  regardless  of  their  speed 
range  or  medium  of  operation.  Without  such  knowledge,  vehicle  per- 
formance cannot  be  adequately  predicted  and  considerable  overdesign 
is  mandatory  to  insure  that  the  veliicle  will  satisfy  the  design  criteria. 
Ship  hulls,  aircraft,  reentry  vehicles,  and  the  space  shuttle  are  but 
a few  areas  in  which  the  design  is  influenced  heavily  by  the  knowledge 
or  lack  of  knowledge  of  roughness  effects. 

One  of  the  earliest  extensive  studies  of  roughness  effects  was 
that  of  Nikuradse,^  in  which  the  effects  of  Reynolds  number  and 
relative  roughness  heights  were  investigated.  This  work  w.is 
carried  out  using  dense  sand  grain  roughness  in  pipe  flow.  He  dis- 
covered that  for  Reynolds  number  based  on  roughness  height  below  a 
certain  critical  value,  the  roughness  had  no  discernible  effect  on  the 
flow,  i.  e.  , the  surface  was  "hydraulically  smooth."  For  roughness 


Reynolds  numbers  above  this  critical  value  and  below  a second  c ritical 
value,  the  effect  of  the  w.all  roughness  on  the  boundary  lay«;r  properties 
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was  dependent  on  flow  Reynolds  number  and  roughness  density,  k/d 
(k  is  the  rouglmess  height,  and  d is  some  characteristic  length  of 
the  flow,  such  as  the  pipe  diameter  for  Nikuradse's  work).  Tlutse 
flows  were  termed  "transitional.”  Flows  with  rovighness  Reynolds 
numbers  above  the  second  critical  value  were  found  to  depend  only 
on  the  roughness  density  and  were  termed  "fully  rough." 

Moore, in  1951,  investigated  a zero  pressure  gradient 
boundary  layer  over  a roughness  consisting  of  square  bars  placed 
normal  to  the  flow  with  a ratio  of  pitch  to  height  of  4.  Moore  found 
that  a similarity  defect  law  correlated  his  boundary  layer  profiles, 
and  the  law  was  identical  with  the  smooth-wall  law,  provided  the  origin 
for  measuring  y and  6 (the  boundary  layer  thickness)  was  located 
some  distance  below  the  crest  of  the  roughness  elements. 

(3 ) 

Hama,'  ' in  1954,  conducted  an  extensive  investigation  which 
(4) 

showed  that  the  Clauser'  ' form  of  the  logarithmic  velocity  distribution 


for  rough  wall  flows 


where 


u 1 ^ T 


= Vt  / p 

• X*/  • It 


1C  and  C 


universal  constants 


wall  shear  stress 


= fluid  density  at  the  wall 


a roughness  function  which  is  zero  for  smooth 

walls  and  which  depends  on  the  roughness  Reynolds 
number 
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and  the  Clauser'  ' form  of  Uie  roughness  function  for  fully  rovigh 


flows 


A 1 h u 

*1!  = >(— !)  I I) 
u K ' V ' 

T 


are  both  universal  for  a given  roug.hness  gef>nu‘lry  in  pipe,  channel^ 
and  zero  pressure  gradient  boundary  layer  flow.  The  constant  1)  is 
found  to  depend  upon  the  free  stream  pressvire  gradient. 

More  recently,  the  work  of  Betterman,^^^  Morris, I.iu 
et  al.  and  Perry  et  have  included  the  effects  on  the  constant 

D of  the  element  density.  Perry  and  Joubert  have  investigated  the 
effect  of  an  adverse  pressure  gradient  on  the  rougliness  function.  All 
of  these  works  were  for  either  a sand  grain  roughne ss^  ^ ^ or  a trans- 
verse square  bar  type  of  roughness, ^ while  Streeter  and  Chu,^^^^ 
Sams,^^^^  Ambrose, and  Corrsin  and  Kistler,^^^^  have  investigated 
different  types  of  roughnesses. 

Antonia  and  Luxton,^^"^^  Liu  et  al.  and  Logan  and  Jones^^^’^ 

are  among  those  who  have  stvidied  the  turbulent  prop«!rtii-s  of  incom- 
pressible fluid  flow  over  rough  surfaces.  They  have  found  that  the 
velocity  fluctuation  magnitvides  in  the  outer  jjart  of  the  rough  wall 


boundciry  layer  are  significantly  higlier  than  those  on  the  smooth  wa’l. 
The  fluctuation  profiles  w'ere  also  found  to  .issume  a self  preserving 
shape. 

Using  the  information  obtained  in  these  investigations, 

(17) 

Van  Driest  has  constructed  a matheniatical  model  of  incompres- 


sible fluid  flow  over  a rough  wall  which  has  met  with  considerable 
success  when  used  to  compute  the  mean  properties  of  such  a flow. 
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Dvorak^^^^  utilizt^s  the  correlation  ol  lletterruan,  extendo  it  to  otlier 
roughness  densities,  and  incorjior ater.  tlie  universal  rough  surface  law 
of  the  wall  into  a proc^;dure  for  computing  incompressible  rough  wall 
turbulent  boundary  layer  profiles. 

Considerably  less  progress  has  been  made  in  determining  the 
effects  of  surface  roughness  on  compressible  turbulent  boundary 

(19) 

layers.  One  of  the  first  inve.stigations  was  that  of  Goddard  who 
studied  the  effects  of  sand  grain  roughness  .at  M.ich  numbers  c)f  2 to 
4.  9.  Wade^^*^^  at  about  the  sanae  time  determined  the  i fleets  of  "screw 
th^e^ld"  roughness  at  Mach  nuniber  2.  S.  Fenter,  Young,  Shutts 
and  Fenter,  Mann,  Red.a,  and  Monta  et  al.  h.ave  .also 

investigated  this  phenomenon  at  Mach  numbers  varying  from  2 to  9. 

With  the  added  phenomenon  of  compressibility,  it  is  no 
longer  possible  to  directly  correlate  the  boundary  layer  profiles 
in  the  law-of-the-wall  form  tis  was  done  for  incompressible  flow. 


Even  compressible  srnootli  wall  turbulent  boundary  Layers,  at  dif- 
ferent Mach  numbers  or  wall-to-stagnation  temperature  r.atios. 


cannot  be  directly  correlated  by  the  law  of  the  wall.  Several  trans- 
(27- 37) 

form  methods'  ' have  been  proposed,  based  on  either  theoretical 


analysis  or  experimental  dald. 


transform  compressible  velocity  profile  data  into  an  "equivalent" 
incompress  ilile  form  which  can  then  be  correlated  by  the  incompress- 
ible law  of  the  wall. 

Spalding  and  Chi,  Hopkins  et  al.  , ^ Miles  and  Kim, 
and  Hopkins  and  Inouye^^^^^  liave  compared  sever.al  of  the  theories 
with  available  d.ita  and  conclude  that  while  some  of  them  work 
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1 


quite  well  for  a wide  variety  of  data,  many  of  them  perform  ade- 
quately only  for  data  which  were  acquired  under  certain  conditions 
such  as  limited  Mach  number  range,  lirnited  wall  temperature 
range,  etc.  Since  only  the  immediate  wall  region  of  the  compress- 
ible turbulent  boundary  layer  is  affected  by  the  presence  of  wall 
roughness,  it  seems  likely  that  the  smooth  wall  transform  methods 

will  work  for  the  rough  wall  case  as  well,  and  Fenter,  Young, 

(2  5) 

and  Reda  have  verified  this. 

Dvoraik^  ^ and  Chen^  have  developed  computational 
procedures  based  on  a combination  of  empirical  laws  and  theo- 
retical equations  which  purport  to  predict  the  compressible  rough 
wall  boundary  layer  development  in  a pressure  gradient.  The 
agreement  with  the  limited  amount  of  data  is  quite  good,  but  Dvorak 
specifically  deplores  the  availability  of  suitable  test  data. 

The  earliest  comprehensive  investigation  of  the  fluctuating 

properties  of  the  compressible  turbulent  boundary  layer  was  per- 
(43) 

formed  by  Kistler'  . He  found  that  in  the  Mach  number  range 
of  1.7  to  4.7,  the  mass  flow  and  total  temperature  fluctuation 
intensities  increased  throughout  the  boundary  laycm  with  increasing 
Mach  number.  The  velocity  fluctuation  profile  was  observed  to  be 
generally  similar  in  shape  to  low  speed  fluctuation  profiles.  Sub- 
sequent to  Kistler' s work,  Owen  and  Horstman*'^'^^  and  Laderman 
and  Demetriades,  have  investigated  smooth  wall  turbulent 

boundary  layer  flows  at  Mach  numbers  of  7,  8,  and  9 with  adiabatic 
and  cold  walls.  Laderman  and  Demetriades*"*^^  conclude  that  the 
fluctuation  intensity  is  strongly  dependent  on  the  wall  temperature 
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to  stagnation  temperature  ratio. 

An  investigation  of  the  response  of  a hypersonic  turbulent 
boundary  layer  to  a step  change  in  roughness  was  seen  as  serving 
several  purposes;  it  wovild  provide  information  on  the  distance  re- 
quired to  achieve  equilibrium  rough  wall  flow  in  hypersonic  flow-- 
essential  for  future  roughness  work;  it  woidd  yield  information  on 
the  equilibrium  rough  wall  boundary  layer  flow  downstream  of  the 
transition  region,  creating  additional  test  data  for  the  existing  com- 
putational procedure;  and  it  would  serve  as  a preliminary  step  for 
possible  later  investigations  of  a turbulence  production  dominated 
region  of  flow,  designed  to  gain  further  understanding  of  the  mecha- 
nisms involved  in  turbulence.  In  addition,  comparison  of  the  data 

( 14) 

with  the  existing  low  speed  results  of  Antonia  and  Luxton  would 
reveal  the  Mach  number  dependence  of  the  boundary  layer  response 
to  sudden  perturbations. 


A detailed  account  of  the  present  investigation  is  presented 
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DESCRIPTION  OF  THE  EXPERIMENT 

Facility 

This  work  was  performed  in  Leg  II  of  the  Graduate  Aeronautical 
Laboratories,  California  Institute  of  Technology  (GALCIT)  hypersonic 
wind  tunnel.  The  tunnel  is  a closed  circuit,  continuously  operating 
facility  utilizing  heated  air  as  the  test  gas.  Leg  II  has  a two  dimen- 
sional, flexible  nozzle  which  may  be  contoured  to  produce  Mach 
numbers  in  the  range  of  6 to  9.  The  side  walls  of  the  tunnel  diverge 
to  compensate  for  boundary  layer  growth  in  the  flow  direction.  For 
this  work,  a half  nozzle  configuration  was  contoured  to  produce  a 
nominal  free  stream  Mach  number  of  6.  0.  This  resulted  in  a total 
test  section  height  of  approximately  2.8  inches,  with  an  inviscid  core 
approximately  0.8  inches  in  height  by  5.0  irches  in  width. 

The  nominal  tunnel  operating  conditions  for  this  work  were 
chosen  to  give  the  highest  possible  Reynolds  number  consistent  with 
good  quality  flow  and  safety  conditions.  These  conditions  were  p^  = 228 
psia  and  T^  = 770°R  which  resulted  in  a free  stream  Reynolds  number. 
Experimental  Apparatus 

The  bottom  flexible  nozzle  wall  of  the  tunnel  was  removed  and 
replaced  by  an  assembly  consisting  of  a permanent  0,  625  inch  thick 
steel  plate  (base  plate)  upon  which  removable  surface  plates  0.  35  inch 
thick  were  mounted.  The  base  plate  itself  was  48.  8 inches  long  and 
the  removable  surface  element  length  was  46.  3 inches. 

Four  sets  of  surface  plates  were  fabricated,  one  of  which  was 


left  smooth  to  serve  as  a reference  condition,  while  the  others  were 
machined  to  different  surface  roughness. 
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The  surface  roughness  chosen  for  this  study  was  a transverse 
square  bar  type  as  illustrated  in  Figure  1.  This  configuration  was 
chosen  because  it  has  been  extensively  studied  in  low  speed  flow 
(Moore,  Antonia  and  Luxton,  Perry,  Schofield  and  Joubcrt, 
Betterman^^^)  and  it  may  be  readily  characterized  by  a height  and  a 
wavelength.  The  wavelength  to  height  ratio  used  was  4,  and  the 
grooves  were  recessed  into  the  surface,  leaving  the  top  of  the  rough- 
ness elements  level  with  the  smooth  wall  upstream  of  the  roughness. 

Roughness  heights  of  0.01Z5,  0.0250  and  0.050  inch  were  used,  which, 

(5) 

using  the  data  of  Betterman,  correspond  to  incompressible  sand 
grain  roughness  heights  of  0.  04,  0.  09,  and  0.  18  inch  yielding  non- 

*1  S T 

dimensional  roughness  values  k (=  — ) of  approximately  18,  40 

w 

and  85. 

To  reduce  the  heat  loss  from  the  flat  plate  to  the  room,  a con- 
tinuous flow  of  low  speed  (<  40  ft/sec)  heated  air  was  maintained  in  a 
channel  formed  by  the  base  plate  and  the  tunnel  side  walls. 
Instrumentation 

The  mean  flow  properties  of  the  boundary  layer  were  computed 
from  data  acquired  during  Pitot  pressure  and  total  temperature  surveys. 
A limited  amount  of  static  pressire  data  was  obtained,  the  main  concern 
being  the  verificiation  of  the  standard  assumption  of  constant  static 
pressure  across  the  boundary  layer.  Direct  measurements  of  wall 
temperature  and  pressure  were  made  via  instrumentation  in  the  tunnel 
wall,  while  the  wall  shear  stress  was  determined  using  a skin  friction 
balance  which  could  be  installed  at  either  of  two  axial  locations.  A 
constant  current  hot  wire  anemometer  was  used  to  obtain  fluctuation 
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data  from  which  the  turbulent  properties  of  the  boundary  layer  could 
be  deduced, 

Pitot  pressures  were  measured  using  a Statham  pressure  trans- 
ducer (PA-208TC- 10-3 50,  0-10  psia).  Static  pressures  were  measured 
using  a Datametric  electronic  barometer,  t'^'pe  1014A.  Pitot  pre s sure 
readings  are  estimated  to  be  accurate  to  within  :r 

The  total  temperature  probe  was  based  on  tlie  design  of 

(47) 

Behrens'  and  consisted  of  a micarta  body  and  wedge  shaped  micarta 
supports  0.  5 inch  apart. 

Two  traverse  mechanisms  were  used  in  this  experiment.  One 
traverse  was  window  mounted  and  permitted  cross  stream  as  well  as 
vertical  and  axial  movement,  while  the  other  was  a top  mounted  unit 
which  utilized  a more  rigid  support  system  for  the  axial  moving  sting 
and  permitted  only  axial  and  vertical  movement.  Both  mechanisms 
were  located  at  the  aft  end  of  the  test  section  (a;  52  inches  from  the 
throat).  Estimates  of  the  probe  position  (with  respect  to  the  flat  plate) 
accuracy  under  tunnel  operating  conditions  vary  from  ± 0.  002  inch  at 
the  aft  end  of  the  test  section  to  i 0,  005  inch  at  the  maximum  probe 
extension  of  28  inches. 

The  smooth  surface  plate  sections  were  instrumented  with 
static  pressure  taps  and  thermocouples.  The  taps  consisted  of  a 0.  067 
inch  diameter  hole  drilled  from  the  back  surface  of  the  surface  plate 
to  within  0.025  inch  of  the  front  surface  with  a 0.014  diameter  inch 
hole  through  the  surface.  These  static  pressure  readings  are  estimated 
to  be  accurate  to  within  ± 1 mm  silicone  oil  (0.0013  psi)  or  ± T^’  or  less 
of  the  minimum  wall  pressure  sensed.  The  outputs  of  all  of  the  surface 
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thermocouples  were  monitored  on  a self  balancing  strip  chart  recorder, 
and  the  accuracy  of  the  thermocouples  based  on  calibration  was  ± 1*^F. 

Provisions  were  made  for  mounting  a skin  friction  balance  at 
two  axial  locations,  27.9  and  47.9  inches  from  the  throat.  The  surface 
elements  of  the  balance  could  be  interchanged  in  the  same  manner  as 
the  rest  of  the  flat  plate  surface  sections.  Roughness  elements  for  the 
balance  were  machined  to  match  the  various  plate  roughnesses  with  no 
interruption  or  discontinuity.  The  roughness  heights  were  such  that  an 
integral  number  of  roughness  wavelengths  would  be  present  on  the  bal- 
ance element  (1,  2,  and  4X  for  k=  0.050,  0.025,  and  0.125  inch,respectively). 

The  balance  itself  was  a floating  element,  null  return  instrument 
based  on  the  design  of  Coles.  Installation  of  the  instrument  in  the 

flat  plate  was  accomplished  by  means  of  lb  leveling  screws  located 
near  the  outer  edge  of  the  balance  table  which  were  alternately  tapped 
into  the  table  and  into  the  base  plate.  Discontinuity  in  surface  level 
between  the  balance  table  surface  and  the  surrounding  plate  surface 
was  held  to  less  than  0.0001  inch  for  the  smooth  flat  plate  and  less  than 
0.  0002  inch  for  the  rough  plates. 

The  hot  wire  anemometer  system  used  in  this  experiment  was  a 
Shapiro- Edwards  constant  current  set  with  a half  power  frequency  of 
320  KHz. 

Experimental  Procedure 

The  nozzle  contour  was  set  with  the  smooth  surface  to  provide 
uniform  Mach  6 flow  in  the  test  section,  and  this  contour  was  not  altered 


when  different  surfaces  were  used.  A nominal  contour  was  calculated 
using  a method  of  characteristics  computer  program  and  correcting 


for  the  boundary  layer  displacement  thickness.  This  contour  was  then 
adjusted  (while  the  tunnel  was  running  at  operating  conditions)  to  mini- 
mize freestream  Mach  number  variations  in  tlie  test  section.  The 
centerline  axial  Mach  distribution  achieved  in  this  manner  is  shov'n 
in  Figure  2. 

The  surface  plate  lengths  were  such  that  the  leading  edge  of 
the  third  plate  from  the  throat  was  located  at  x = 25.  4 - several 
boundary  thicknesses  downstream  of  the  intersection  of  the  test 
rhombus  with  the  wall.  All  step  changes  in  surface  roughness  occurred 
at  this  location.  For  the  smooth  surface  to  rough  surface  configurations, 
the  surface  of  the  two  aft  surface  plates  consisted  of  the  roughness 
illustrated  in  Figure  1 with  one  of  three  roughness  heights  (k  = 0.0125, 
0.025  and  0.050  inch).  Preliminary  results  of  the  smooth-to- rough 
wall  investigation  revealed  that  an  axial  distance  of  some  8 inches 
(106  or  206^)  w’as  sufficient  for  the  disturbed  boundary  layer  to  reach 
a new  mean  flow  equilibrium  state  over  the  rough  wall.  Based  on  this 
information,  the  rough-to- smooth  wall  step  change  model  was  con- 
figured with  a smooth  plate  immediately  downstream  of  the  throat, 
followed  by  rough  surface  (k  = 0.050  inch)  second  plate,  with  the  use 
of  two  smooth  surface  plates  downstream  of  x = 25.  4 inches  completing 
the  surface  plate  set.  The  length  of  the  rough  wall  section  was  11.5 
inches  - considerably  more  than  the  length  found  necessary  for  the 
attainment  of  new  mean  flow  equilibrium  profiles  in  the  absence  of 
an  axial  pressure  gradient.  The  choice  of  k = 0.050  inch  for  the 
roughness  height  was  predicated  on  the  fact  that  preliminary  calcula- 
tions indicated  it  would  be  in  the  fully  rough  regime,  while  the  other 
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roughness  heights  would  fall  into  the  transitionally  rough  region. 

The  low  speed  channel  below  the  flat  plate  was  instrumented 
with  two  thermocouples  on  the  lower  side  of  the  base  plates  on 
the  plate  centerline,  two  on  the  lower  side  of  the  base  plate 
adjacent  to  a side  wall,  and  two  in  the  low  speed  flow,  two  inches 
below  the  base  plate.  With  the  tunnel  at  operating  conditions  the 
low  speed  heater  temperature  was  adjusted  so  that  the  tempera- 
tures indicated  by  the  two  thermocouples  on  the  low  speed  side 
of  the  base  plate  were  the  same  as  those  indicated  by  their 
respective  surface  thermocouples.  The  thermocouple  readings 
were  monitored  over  a period  of  several  hours  and  were  found 
to  remain  constant. 

The  thermocouples  near  tlie  side  walls  indicated  tempera- 
tures within  2°F  of  those  on  the  tunnel  centerline,  indicating 
negligible  heat  loss  to  the  side  walls.  There  was  a somewliat 
more  severe  axial  temperature  gradient  near  the  aft  of  the  test 
section.  This  gradient  was  due  to  the  presence  of  the  support 
structure  for  the  aft  end  of  the  flat  plate  and  could  not  be  elim- 
inated. The  influence  of  this  slight  gradient  on  the  plate  tem- 
perature is  considered  to  have  been  negligible  as  evidenced  by 
the  very  long  section  of  constant  plate  temperature  shown  in 
Figure  3. 

A simple  calculation  based  on  the  tunnel  stagnation  tem- 
perature indicates  that  a recovery  temperature  of  695°R  would 
be  expected.  As  shown  in  Figure  3,  the  flat  plate  temperature 
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was  618°R,  considerably  below  the  anticipated  value.  The  reason 
for  this  discrepancy  is  not  fully  understood,  but  it  is  believed  to  be 
due  to  a conibination  of  heat  loss  from  the  uninsulated  upper  nozzle 
wall,  and  the  influence  of  throat  cooline;  far  upstream.  Tliis  will  be 
discussed  further  in  the  Results  and  Discussion  section,  below. 

Prior  to  the  acquisition  of  any  data  the  flat  plate  was  brought 
up  to  within  2-3°F  of  its  equilibrium  temperature  over  its  entire 
length.  The  temperature  distribution  on  the  plate  surface  v,ould 
normally  level  out  and  reach  equilibrium  conditions  some  45  minutes 
to  1 hour  after  tunnel  startup. 
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DATA  REDUCTION 

Pressure  Data 

The  free  stream  Mach  number  was  computed  from  the  free 
streani  total  pressure  to  stagnation  pressure  ratio  assuniing  an 
isentropic  expansion  from  the  throat  to  the  final  flow  conditions. 

The  static  pressure  corresponding  to  this  Mach  number  and  stag- 
nation pressure  was  then  computed  (using  isentropic,  perfect  gas 
relations)  and  assumed  to  be  constant  across  the  boundary  layer 
(a  static  pressure  traverse  established  ttie  validity  of  this  approxi- 
mation). The  local  measured  Pitot  pressure,  corrected  for  Reynolds 
number  effects  using  the  results  of  Ramaswamy , ^ was  used  in 

conjunction  with  tins  static  pressure  to  compute  the  local  Mach 
number  from  the  Rayleigh  Pitot  formula. 

Total  Temperature  Dat<. 

The  local  total  temperature  was  determined  using  the  theory 
of  Behrens. Due  to  the  probe  construction,  it  was  not  possible 
to  acquire  data  within  0.  030  inch  of  the.  wall,  so  a linear  variation 
of  total  temperature  between  the  last  measured  data  point  and  the 
wall  temperature  was  assumed.  Adiabatic,  perfect  gas  relations 
were  used  to  compute  the  local  static  temperature. 

Hot  Wire  Data 

Quantitative  information  on  the  fluctuating  flow  properties 
was  obtained  from  the  hot  wire  measurements  utilizing  the  techniques 
developed  by  Kovasznay , ^ ^ Morkovin,^^^^  Kistler,^  ^ Laufer,^  ^ 

and  Gran.*"^*^^  The  reduction  of  both  mean  and  fluctuating  flow  hot 

wire  data,  including  the  end  loss  correction,  is  given  in  a condensed 
(49) 

form  by  Gran. 
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RESULTS  AND  DISCUSSION 


Profile  Data 

Mean  flow  Pitot  pressure  and  total  temperature  data  were 
obtained  throughout  the  boundary  layer  at  least  every  two  inches 
axially  throughout  the  test  section.  Near  llie  step  change  in  roughness, 
profile  data  were  obtained  at  one  inch  axial  spacings. 

Smooth  Wall 

Comparison  of  the  smooth  wall  velocity  profiles  as  presented 
in  Figure  4 reveals  that  the  smooth  plate  boundary  layer  is  very 
nearly  self  preserving  downstream  of  x = 25.4.  A small  amount  of 
adjustment  in  the  shape  is  seen  to  occur  between  x = 2 5.  4 and  x = 29.  4. 
followed  by  very  nearly  identical  profiles  downstream  to  x = 3 5.  4. 
However,  the  velocity  profile  suddenly  becomes  more  "full''  near  the 
wall  at  X = 27.  4.  Although  no  data  downstream  of  x = 37.  4 are  pre- 
sented on  the  overlay  plot,  the  profile  data  remain  self- similar  down- 
stream to  X = 47.4.  The  velocity  profiles  for  all  the  step  change  con- 
figurations exhibit  similar  behavior  in  the  vicinity  of  x = 37.  4. 
Comparable  perturbations  are  present  in  the  density  data.  The  profile 
integral  data  (to  be  discussed  later)  in  Figures  8,  9,  10  and  11  also 
exhibit  sudden  shifts  in  level  in  the  vicinity  of  x = 37.4. 

An  examination  of  the  typical  freestream  Pitot  pressure  distri- 
bution (Fig.  2)  offers  some  insight  into  the  cause  of  these  anomalies. 
Small  pressure  peaks  are  seen  to  be  centered  in  the  freestream 
(y  = 1.1  inches)  near  x = 33  and  x = 48.  This  spacing  is  compatible 
with  the  existence  of  a weak  pressure  wave  which  is  reflected  off  the 
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upper  nozzle  wall  such  that  it  intersects  the  lower  wall  bounclary 
layer  near  x = 33,  reaches  the  wall  near  x = 40  iind  reflects  back 
into  the  freestream  near  x = 47.  Static  pressure  Iroiti  widely  dis- 
persed smooth  plate  boundary  layer  surveys  indicate  that  such  a 
pressure  wave  does  exist,  although  it  is  weak -- maximum  pressure 
variation  across  the  boundary  layer  was  observed  to  be  some  1^ 


at  axial  location  3 3.4. 

Since  all  the  data  were  reduced  assuming  = 0,  a question 
arises  as  to  whether  the  observed  effects  cif  the  pressure  wave 


are  real  or  due  to  the  method  of  data  redxiction.  Reduction 


of  tlie  smooth  plate  profile  data  at  x = 31.4  and  39.4  utilizing 
the  measured  static  pressure  distributions  produced  velocity  and 
density  profiles  which  were  very  similar  between  the  two  locations. 
The  effect  of  the  wave  was  observed  to  be  a slightly  fuller 
velocity  profile  close  to  the  wall  at  x = 33.  4 as  compared  to 
X = 39.  4.  The  effect  of  this  wave  on  the  profile  integral  data 
will  be  discussed  in  a later  section.  Thus,  even  though  the 
pressure  variation  across  the  boundary  layer  was  small,  the 
assumption  = 0 is  the  major  cause  of  the  non  self- similar 
boundary  layer  profiles.  The  presence  of  the  pressure  wave 
does  affect  the  profiles,  but  only  slightly.  Use  of  the  measured 
static  pressure  data  was  precluded  by  the  limited  amount 
obtained  (some  four  axial  locations,  all  on  the  smooth  plate), 
so  the  constant  pressure  assumption  was  utilized  to  produce  the 
data  reported  here. 
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Comparison  of  the  profile  data  for  the  various  configurations 
was  predicated  on  the  fact  that  although  strictly  similar  profiles 
could  not  be  attained  for  an  extended  axial  distance  due  to  the 
pressure  wave,  the  effect  of  this  wave  on  tlic  profile  should  be 
independent  of  configuration. 

Roughness  Effects 

The  step  change  configuration  profih'  data  are  presented 
in  Figures  5,  6 and  7.  A cursory  examination  of  the  profile 

data  reveals  that  tlie  step  change  in  rouglmess  does,  as  expected, 
introduce  significant  changes  into  the  various  profiles  and  results 
in  new  equilibrium  or  self- similar  profiles  some  distance  down- 
stream. 

The  development  of  the  disturbed  boundary  layer  may 
readily  be  observed  from  overlays  of  profile  data  at  progressively 
greater  distance  downstream  of  the  throat  (Figs.  5,  6 and  7). 

The  smootli-to- rough  wall  development  as  seen  in  these  figures 
is  representative  of  all  the  smooth-to- rough  configurations.  An 
increase  or  decrease  in  roughness  size  simply  causes  a corre- 
sponding increase  or  decrease  in  the  magnitude  of  the  observed 
effects.  The  influence  of  the  step  change  in  roughness  is  seen 
to  spread  rapidly  across  the  entire  layer,  causing  progressively 
larger  changes  as  the  layer  continues  downstream,  until  a new 
equilibrium  profile  is  assumed. 
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Analysis  of  the  srnooth-to- rough  configuration  profiles  reveals 

that,  for  all  quantities  presented,  the  attainment  of  a new  equilibrivim 

profile  is  accomplished  in  approximately  the  same  axial  distance:, 

regardless  of  the  roughness  height.  Tlie  neu  profiles  are  definitely 

established  inches  downstream  of  the  throat,  cor respcjnding  to 

some  10  boundary  layer  thicknesses  (6)  or  ZO  displacement  thick- 
* 


nesses  (6  ) downstream  of  the  step  change. 

The  boundary  layer  appears  to  take  a slightly  greater  distance 

to  adjust  to  the  rough-to- smooth  step  change,  but  even  in  this  case 

the  new  equilibrium  profile  is  definitely  established  37  inches  down- 

stream  of  the  throat,  some  146  or  Z66  downstream  of  the  step 

change.  Thus,  although  the  boundary  layer  adjusts  to  the  rough- to- 

smooth  step  change  somewhat  more  slowly  than  to  tlie  srnooth-to- rough 

step  change,  the  difference  in  distance  involved  is  only  some  30,'^  in 
* . 

terms  of  6 , or  in  terms  of  6. 

t 54) 

Jacobs  performed  his  low  speed  work  in  fully  developed 
channel  flow  on  a two-dimensional  roughness  comprised  of  transverse 
rectangular  bars.  He  found  that  the  rough-to- smooth  surface  change 
resulted  in  a new  eejuilibrium  state  within  17  channel  half  heights, 
while  the  rough- to- smooth  change  took  some  25  channel  half  heights, 
an  increase  of  50)b.  Antonia  and  I.uxton,  ' working  with  the 
same  type  of  roughness  as  used  liere  in  a low-speed  boundary  layer, 
found  that  less  than  ?.0  boundary  layer  thicknesses  were  required 
for  the  flow  over  the  smooth-to- rough  step  to  .idapt  to  the  change, 
while  after  16  boundary  layer  thicknesses  the  profiles  were  "far  from 
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self  preserving"  for  the  flow  ovtn-  Ihe  rougli- lo- smooth  step. 

The  present  data,  when  6 is  regarded  as  tiie  tliickness  of  the  com- 
pressible boundary  layer,  are  consistent  with  the  results  of  both 
investigations. 

Comparison  of  the  profile  data  for  the  smooth  wall  case  and 
the  rough-to- smooth  step  change  c.ise  (as  jiresented  in  Figure  7) 
reveals  that  although  the  flow  over  the  step  change  does  attain  new 
equilibrium  profileis,  they  are  not  the  same  as  the  profiles  for  the 
smooth  wall.  The  velocity',  density  and  Mach  number  profiles  for 
the  rough-to- smooth  configuration  are  found  to  be  considerably  less 
full  throughout  the  boundary  layer  than  the  corresponding  smooth 
plate  configuration  profiles.  Much  smaller  differences  are  noticeable 
in  the  total  temperature  data.  Thus  the  presence  of  the  rough  wall 
seems  to  give  rise  to  some  non- reversible  ch.mges  in  the  boundary 
layer  flow. 

Rotta^'’^^  and  Bertram  and  Neal/^^'  have  suggested  that  non- 
equilibrium or  upstream  history  effects  may  account  for  the  differ- 
ences observed  in  nozzle  wall  and  flat  phite  boundary  layers. 

Feller , ' Voisinet,  et  al.|^^^  Stur<-k^^*^^  and  Bushnell  , 
et  al.^^^^  have  investigated  the  effects  of  upstream  conditions  on  the 
bmmdary  layer.  Their  findings  indicate  that  the  properties  of  the 
nozzle-wall  boundary  layers  are  indeed  very  sensitive  to  changes 
in  the  upstre.irn  wall  temperature  and  freestream  pressure  gradients. 
The  effect  of  the  upstream  favorable  pressure  gradient  was  found  by 
Bushnell  , ot  al.  to  result  in  increased  fullness  of  the  velocity  profile 
while  not  significantly  affecting  the  total  temperature  profile.  They 
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quote  the  data  of  Feller  and  Jones  in  the  Langley  Mach  b liigh  Reynolds 
number  tunnel  as  revealing  a tendency  to  relax  from  the  characteristic- 
tunnel  wall  quadratic  total  temperature-velocity  variation  to  the  usual 
flat  plate  linear  relationship  some  60  boundary  layer  thicknesses 
do-wnstream  of  the  nozzle  exit.  If  the  data  of  Feller  and  Jones  could 
be  applied  to  this  work,  this  relaxation  tendency  would  become  evident 
some  three  feet  downstream  of  the  end  of  the  test  section. 

If  the  subjection  of  the  boundary  layer  to  the  rough  wall  section 
followed  by  the  readjustment  to  smooth  wall  conditions  produced  a large 
enough  perturbation  to  destroy  or  significantly  decrease  these  upstream 
pressure  and  temperature  effects,  the  anticipated  result  would  be  a 
somewhat  less  full  velocity  profile  and  a more  nearly  Crocco  Temperature- 
velocity  variation.  The  observed  velocity  profile  difference  supports 
such  a hypothesis,  as  do  the  temperature-velocity  relationships  pre- 
sented in  Figure  IZ. 


Profile  Integral  Data 
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Jhese  data  are  tabulated  in  Table  11  and  presented  graphically  in 
Figures  8,  9,  10  and  11. 
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As  discussed  earlier,  the  flow  in  the  test  section  where  these 
measurements  were  obtained  was  characterized  by  a freestream 
Mach  number  of  6.  02  ± 1^'  where  the  variations  in  the  Mach  number 
were  found  to  be  caused  by  tlie  presence  of  a weak  pressure  wave. 

The  effect  of  this  wave  on  the  profile  data  was  discussed  earlier  with 
respect  to  the  consequences  of  utilizing  the  assumption  = 0 for 

the  boundary  layer  reduction.  The  effect  of  tiiis  assumption  on  the 
integral  data  is  illustrated  in  Figures  8,  9,  10  and  11  where  the  sum- 
mary parameters  resulting  from  the  use  of  the  measured  static  pres- 
sure profile  at  two  axial  locations  on  the  smooth  wall  are  plotted. 

The  = 0 assumption  is  observed  to  be  the  major  reason  for 
the  anomalous  behavior  of  the  momentum  thickness  near  x = 35.  This 
assumption,  however,  has  little  effect  on  the  displacement  thickness 
and  enthalpy  flux,  implying  that  their  sudden  shift  in  level  in  the  vicinity 
of  X = 33  — 37  is  due  to  the  presence  of  the  weak  pressure  wave  in  the 
test  section. 

The  momentum  integral  equation  for  a two  dimensional,  zero 
pressure  gradient  flow  such  as  this  may  be  written 


S 

dx  ~ 2 ~ 2 

"^e  e 

The  skin  friction  was  directly  measured  (through  use  of  the  sk-n  friction 
balance)  for  three  configurations  (smooth  plate  and  smooth-to- rough 
wall  with  k = 0.  025  and  0.  050  inch)  at  x = 27.  9 and  47.  9.  These  data 
are  included  in  Table  II.  The  data  at  x = 47.9  were  used  in  the  above 
equation  to  determine  the  anticipated  streamwise  variation  of  0 for  the 
various  configurations.  Lines  with  the  appropriate  slopes  are  faired 
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through  the  data  downstream  of  x = 3 7 (where  the  presence  of  the  pres- 
sure wave  has  little  effect  on  the  data)  in  Figure  9.  The  actual  stream- 
wise  variation  of?  is  seen  to  be  in  good  agreement  with  the  calculated 
value  for  all  configurations. 

Examination  of  the  data  in  Figures  8-11  with  the  effect  of  the 
wave  on  the  data  in  the  vicinity  of  x = 3 3-37  in  mind  reveals  no  sudden 
changes  due  to  the  existence  of  a step  change  in  roughness.  The 
parameters  are  observed  to  adjust  gradually  to  the  change  and  tcj 
approach  new  equilibrium  behavior  in  the  vicinity  of  x = 37,  just  as 
the  mean  flow  profiles  do. 


Transformation  and  Correlation 

As  mentioned  in  the  introduction,  there  exists  a multitude  of 
transformation  techniques  for  casting  compressible  velocity  profile 
data  into  an  equivalent  incompressible  form,  but  several  recent 
papers  including  Lewis,  et  al.  , Keener  and  Hopkins, and  Kemp 

and  Owen^^^^  have  found  that  the  Van  Driest  method  accomplishes 

this  task  as  well  as  any  of  the  others,  and  perhaps  better  than  most. 
This  occurs  in  spite  of  the  fact  that  it  was  one  of  the  earliest  developed, 
is  based  on  a less  firm  "theorfetical"  foundation,  and  is  simpler  to  use 
than  most. 

Briefly,  the  Van  Driest  1 transformation  starts  with  the  assump- 
tion of  a Crocco  temperature  distribution  through  the  boundary  layer. 


i.  e. 
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and  uses  the  Prandtl  mixing  length  theory  to  determine  the  Reynolds 
stress  in  the  flow.  The  result  is  a transformation  of  the  compressible 
velocity,  u,  into  an  '’incompressible''  velocity  u''  via  the  eo,uation 


U 

u 

e 


1 ,.i2A(u/uJ-B  ^ 


-r  Sin 

A 


(B  +4A 


+ -r  Sin 

A 


I B 

2 2 ' 
(B  +4A  )‘^ 


(3) 


where  A and  B arc  parameters  that  depend  on  the  temperature  distri- 
bution. 

Van  Driest  deduced  that  this  velocity  would  be  correlated  by 


1 V u 

J_,  3 

< V 


+ F 


where 


1/T  /p 
w ^W 


< = 0.  4 - Karman's  constant 

= kinematic  viscosity  at  the  wall  temperature 
F = constant  dependent  only  on  the  boundary  conditions 

(i.  e.  freestream  pressure  gradient  and  wall 
conditions ) 

The  similarity  of  this  equation  and  the  law  of  the  wall  immediately 
suggests  the  use  of  the  law  of  the  wall  to  correlate  the  transformed 
velocities.  The  form  of  the  law  of  the  wall  to  be  used  in  this  work  is^^^^ 


y'*’  + C + -^  w(y/6)  - 

T 


(4) 


where 
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w(y/6)  = 


Au 
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w 


5.  0 - universal  constant 

wake  strength  parameter  dependent  on  the  frec- 
stream  pressure  gradient 

2 sin^  ” "wake  f\mction"  suggested  by  Coles^^^^ 


Clauser'  s 
k u 


rough  wall  velocity  defect 
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D = a constant  dependent  only  on  the  wall  conditions 

As  stated  above,  the  Van  Driest  transformation  is  based  on  the 
assumption  of  a Crocco  temperature  profile.  The  temperature  profile 
most  often  used  to  estimate  temperature  distributions  in  boundary 
layers  are  one  due  to  Crocco  and  another  given  in  VValz.^^^^  Com- 
parison of  these  formulae  with  the  experimental  data  for  two  sample 
cases  is  made  in  Figure  12  where  it  is  evident  that  while  the  trend 
of  the  temperature-velocity  relation  is  represented  by  both  profiles, 
neither  is  a good  approximation.  Also  shown  is  a least  squares 
quadratic  curve  fit  to  the  data  for  each  case. 

Since  the  Crocco  distribution  is  such  a poor  fit  to  the  data, 


doubt  is  cast  on  the  validity  of  the  transformation  for  the  velocity 
data.  It  is  readily  verified,  however,  that  the  Van  Driest  transfor- 
mation may  be  altered  to  utilize  the  coefficients  of  the  curve  fit  to 
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‘ the  actual  data  in  evaluating  A and  B in  equation  (i).  Both  the  Crocco 

and  empirical  temperature  distributions  were  used  in  an  attempt  to 
correlate  the  data.  The  transformation  of  the  velocity  data  via  the 
"modified"  van  Driest  relationship  resulted  in  a better  correlation 
(in  the  sense  that  tlie  deduced  skin  friction  was  closer  to  the  measured 
value)  for  the  smooth  plate  data  at  x = 47.  This  procedure  yielded  a 
skin  friction  within  I'j^  of  the  measured  value  while  the  use  of  the 
standard  Van  Driest  transformation  resulted  in  a skin  friction  value 
some  lower.  The  modified  form  of  the  transformation  has  been 
used  exclusively  in  tlie  results  discussed  below. 

Figure  13  presents  the  results  of  two  correlations,  both  at 
X = 47.  4 inches  from  the  throat.  The  correlation  of  the  smooth  wall 
data  was  computed,  as  indicated  earlier,  assuming  Au/u.^  = 0 while 
the  rough  wall  correlation  was  based  on  the  measured  skin  friction  at 
that  location.  A comparison  of  Coles'  wake  function  with  the  actual 
data  is  made  in  Figure  14.  In  both  figures  the  fit  of  the  data  is  seen 
to  be  excellent. 

Correlation  Results 

The  various  parameters  resulting  from  the  curvefits  are 
presented  in  Table  II  with  graphical  presententation  made  in  Figures 
15,  16  and  17.  A cursory  examination  of  these  figures  reveals  the 
existence  of  sudden  shifts  in  the  parameter  values  at  axial  locations 
of  33  to  39  inches  from  the  throat.  This  is  the  same  region  in  which 
changes  occurred  in  the  profile  shapes  and  the  integral  properties, 
indicating  that  these  anomalies  are  also  due  to  the  previously 
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discussed  pressure  wave  and  data  reduction  techniques. 


The  wake  strength  parameter  tt  (Fig.  15)  is  observed 
to  vary  betessn  1.1  and  1.6  fpr  all  configurations.  This  is 
significantly  larger  than  the  value  of  0,61  advanced  by  Colcs^^^^^ 
as  being  the  representative  value  for  the  incompressible,  smooth, 
flat  plate  equilibrium  boundary  layer.  Once  again,  the  effect 
of  the  = 0 assumption  in  data  reduction  was  investigated  and 
found  to  be  the  cause  of  the  large  change  in  level  between  x = .53 
and  X = 39  (Fig,  15).  As  stated  above,  only  the  wake  strength 
parameters  for  the  smooth  wall  case  and  those  at  x = 25.  4 and 
47.  4 for  the  step  change  configurations  were  computed  from  the 
curvefit.  All  others  were  estimated  from  tliese  values,  as  dis- 
cussed in  the  Transformation  and  Correlation  section. 


The  existence  of  such  a large  it  parameter  is  in  contrast 
to  the  results  of  Keener  and  Hopkins^^^^  at  Mach  7 (tt  = 0.4) 

and  those  of  Lewis,  et  al.  ' ^ at  Mach  4 (it  = 0,  6).  The  data 

/ / \ 

of  Reda,  et  al.  ' ' at  Mach  2.9  yields  tt  = 1.2,  while  those 

(44)  — 

of  Owen,  et  al.  at  Mach  7 yield  tt  = 0.8.  Laderman  and 

(45)  — 

Demetriades  at  Mach  9 report  tt  = 1.4.  In  all  cases  the 
Van  Driest  transformation  was  used  to  cast  the  velocity  data 
into  an  ''incompressible"  form  for  use  in  the  correlation.  Only 
Keener  and  Hopkins  have  investigated  the  effect  on  the  transfor- 
mation of  using  a fit  of  the  actual  data  in  place  of  the  usual 
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Crocco  distribution.  They  concluded,  as  did  the  author  above, 
that  the  use  of  the  measured  temperature  velocity  relationship 
resulted  in  better  correlation  results  and  a somewhat  higher 
value  for  ir,  compared  to  the  results  adiieved  utilizing  the  Crocco 
distribution. 

The  correlation  skin  friction  results  are  presented  in 
Figure  16,  along  witli  the  measured  skin  friction.  The  measured 
skin  friction  at  x = 47.9  for  the  two  smooth-to- rough  wall  con- 
figurations was  used  to  determine  tr  in  the  correlation,  so 
any  disagreement  of  the  computed  and  measured  values  would 
be  indicative  of  problems  in  the  correlation  procedure.  The 
skin  friction  measurements  at  x = 2 7.9  for  the  step  change 
configuration  and  at  x = 47.9  for  the  smooth  plate,  however, 
were  not  used  as  inputs  to  the  correlation.  The  excellent  agree- 
ment of  the  computed  skin  friction  with  the  measured  values  for 
the  smooth  plate  and  k = 0.025  inch  step-change  configurations 
is  conclusive  evidence  that  the  correlation  procedure  is  valid 
and  that  it  may  be  used  in  a non-equilibrium  region  with  great 
success.  Disagreement  of  the  computed  skin  friction  witli  that 
measured  at  x = 27.9  for  the  k = 0.050  inch  smooth-to- rough 
step  configuration  and  the  behavior  of  the  computed  skin  friction 
in  the  vicinity  of  the  step  indicate  that  the  perturbation  of  the 
Qow  due  to  the  large  roughness  height  is  too  large  for  the 
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correlation  to  be  successful.  The  rough-to- smooth  results^ 
on  the  other  hand,  indicate  a very  smooth  and  relatively  rapid 
transition  to  the  dov-mstream  values.  Although  no  measured 
values  are  available  for  comparison,  th'-st  trends  indicate 
that  the  correlation  is  applicable  for  this  configuration,  oven 
though  the  roughness  height  is  the  same  as  that  in  the  smooth- 
to- rough  step  where  tlie  correlation  was  found  to  be  invalid. 

The  roughness  induced  velocity  shift  data  in  Figure  17 

\ 

(—  is  from  equation  (3)  above)  tend  to  support  tliese  conclusions, 

u, 

although  the  scatter  in  the  k = 0.01Z5  inch  smooth-to- rough 

configuration  would,  by  itself,  give  rise  to  doubts  about  the 

applicability  of  the  correlation  in  this  case.  In  the  light  of 

its  suitability  for  the  k = 0.025  inch  configuration,  however,  this 

scatter  is  attributed  to  the  use  of  the  same  it  (x)  distribution 

for  k = 0.0125  inch  as  was  determined  for  the  k = 0.025  inch 

Au 

configuration.  The  significant  increase  in  the  value  of 
immediately  dowoistream  of  the  step  for  the  k = 0.  050  inch 
smooth-to-rough  case  corresponds  to  either  a significant  decrease 
in  (as  seen  in  Fig.  16)  or  a significant  increase  in  the 
velocity  across  the  boundary  layer.  Neither  of  these  effects 
would  be  expected  and  neither  is  observed  for  the  other  con- 
figurations, thus  substantiating  doubts  about  the  correlation's 
validity  in  the  step  region  for  this  rougliness.  The  rough-to- 
smooth  velocity  defect,  however,  is  observed  to  adjust  gradually 


-29- 


to  the  smooth  wall  conditions,  pausing  at  a somewhat  low  value 
in  the  vicinity  of  x = 37-41  before  attaining  its  smooth  wall 
value  of  0. 

The  examination  of  the  curvefit  results  to  this  point  has 
thus  yielded  the  following  results: 

1.  The  assumed  form  of  the  law-of-the-wall  correlates 
the  equilibrium  smooth  data  very  well,  judging  by  the  accuracy 
with  which  the  skin  friction  was  computed. 

2.  The  correlation  is  applicable  in  the  non-equilibrium 
flow  region  immediately  downstream  of  a step  change  in  rough- 
ness, provided  the  change  is  not  too  large  (what  constitutes 
"large"  cannot  be  deduced  on  the  basis  of  these  results). 

3.  The  boundary  layer  adjusts  gradually  to  the  step 
change,  assuming  the  properties  appropriate  to  the  new  wall 
condition  only  a considerable  distance  downstream  (the  same 
location  as  that  at  which  new  equilibrium  profiles  are  observed). 
This  is  in  contrast  to  the  finding  of  Antonia  and  Luxton^^"^^ 

(for  a smooth-to-rough  step  change)  that  adjustment  to  the 
rough  wall  condition  is  made  within  three  or  four  boundary 
layer  thicknesses  following  an  initial  overshoot  to  a level  above 
that  finally  attained.  Antonia  and  Luxton's^^^^  results  for  a 
rough-to- smooth  change  are  consistent  with  those  found  here, 
i.  e.  changes  gradually  to  the  value  appriate  for  the  smooth 
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Qiiantii:ativf  Rou^liness  Eff ects 

One  convenient  method  of  describing  any  given  ronghne.sa  is  to 
determine  its  "equivalent  sand  grain"  roughness,  which  is  defined  as 
the  height  of  Nikuradse's  sand  grain  roug.hness  which  would  be  required 
to  produce  the  same  velocity  delect.  The  svvhject  rouglmess,  however, 
must  be  "fully  rough,"  i.e.  it  imist  have  a large  enough  roughness 
Reynolds  number  that  the  law  of  the  wall  may  be  expressed  in  the 
form 

u^  = -^2ny/k3  sin^ 

(41 ) 

Dvorak  presents  tentative  criteria  for  the  upper  limit  of  the  tran- 
sitionally smootli  regime  which  indicate  that  tlie  critical  roughness 

, ku 
f T 

Reynolds  number  fo  • tliis  particular  roughness  is  k (=  — — ) - 30. 

•1 

Thus  the  k = 0.  050  inch  roughness  (k  = 33.  8 from  measurements) 
would  be  considered  "fully  rough,"  while  the  k = 0.  0125  inch  and 
k - 0.025  inch  roughness  heights  (k^  = 7.  1 and  14.9,  respectively) 
would  fall  into  the  transitionally  rough  regime. 

For  tile  fully  rough  regime,  the  results  of  Nikuradse  for  in- 
compressible pipe  flo%v  (where  there  is  no  "wake"  function)  may 

i .11  (1,^8) 

be  i-epresented  by  ' 

u = 17  y/k  t 8.5 

K s 

= 

where  is  the  sand  grain  roughness  height.  The  smootli  wall  r«?sults 
of  Nikur.ulse  are  correlated  by 


1 , + 


t 5.  5 


-31- 


Equation  (6)  may  bo  rowr iltoji  a: 
I 


whe  re 


u_ 


- toy  ' 
K 


5.  5 - ^ 
u 

T 


k *■  - 3.0 
K s 


( 7 ) 


vvliich  yields 


K(3  t Ah  ) 
e 


( 8 


As  before,  the  eff(»ct  of  compressibility  is  to  require  the  usf;  of  in 
place  of  V. 

The  application  of  ccjuation  ( 8)  to  the  velocity  defect  of  the 

k=  U.050  inch  roughness  from  Figurt-  33  (AH-  b.ri)yields  k * = 4-1.7  <ir, 

, Ut  s ' 

since  k =33.8,  k /K  = 1.  3.  Returning  for  a moment  to  the  low  speed 


correlation  results  of  Dvor<ik  and  Betterman  (and  extended  by  Siinpson 
and  Dirling^^®^  to  other  geometries  with  goorl  results),  the  effect  of 
the  square  bar  roughne.-;s  with  £/k  = X (waveli-ngth-to-height  ratio) 


(69) 


< 4.  7 is  found  to  he 

Ah 

u 


to  k ' ■(  1 7.  3 S (0.  705  to  X-  1 ) 


(9) 


which  yiehls,  for  k = 33.  8,  X = 4 

Ah  . 8. 4 

u 


The  validity  of  this  equation  is  limited  to  the  fully  rough  regime 
-(  -t 

(k^  70);  howeve  r , k^  =45  is  close  enough  to  this  regime  that  no  large 


errors  are  introduced  ( varies  little  from  the  fully  rough  relation-- 

see  Fig.  18,  This  is  significantly  higher  than  the  observed  value 
(Fig.  1 7 ) of  ().  5-- some  30^  higher.  The  use  of  equation  (8)  with 
•Ah  = 8,  4 yields  k /k  = 4.  8- -more  than  double  the  actual  value  of  1.  3 
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calculated  above.  These  results  illustrate  quite  graphically  the 
inability  of  the  low  speed  correlation  results  of  Dvorak,  Bette rman, 
Simpson  and  Dirling  to  correlate  the  present  liypersonic  data. 

The  current  data  are  presented,  along  with  previously  pub- 
lished data,  in  Figures  18  and  19,  Use  of  the  equivalent  sand  grain 
roughness  Reynolds  number  in  the  abscissa  effects  the  collapse  of 
data  for  several  roughness  configurations  and  Mach  numbers  ranging 
from  0 to  6 into  a single,  well  defined  curve  in  both  cases.  The  sand 

grain  rougliness.  unless  otherwise  noted,  was  computed  for  each  set 

k u,,. 

of  data  by  choosing  a value  of  ^ in  the  fully  rough  regime  and 

w 

utilizing  the  associated  value  of  in  equation  8 above.  This 

value  of  was  then  used  to  plot  all  data  points  for  that  configuration. 

\VTien  presented  in  this  manner,  a single  curve  is  sufficient  to 
describe  the  data  in  Figure  18.  For  k^  > 70  the  appropriate  equation 
is  ti,’!  k^  - 3.  0.  An  equation  of  the  form  suggested  by  Dvorak 

may  be  used  to  describe  the  trend  of  the  data  in  the  transitionally 
rough  regime.  Such  an  equation  is  presented  in  Figure  18.  The 
velocity  defect  and  its  first  derivative  were  specified  to  vanish  at 
k^  = 1.  S and  to  match  the  fully  rough  values  at  k*  = 70, 

The  data  in  Figure  19  are  adequately  represented  by  the  equation 

^f  + 

— ^ = 0.  39(^k  - 2.  3)  + 1. 0 

g.,  s 

'o 

Thus,  once  the  effective  sand  grain  roughness  is  known,  the  effect  of 
surface  roughness  cn  the  skin  friction  and  velocity  defect  may  readily 
be  determined.  The  only  problem  that  remains  is  the  determination 
of  the  equivalent  roughness,  which  is  basically  what  is  accomplished 
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by  tlae  correlations  of  Dirling  and  the  others  mentioned  previously. 
However,  these  correlations  were  found  to  be  inapplicable  for  these 
flow  conditions.  A brief  discussion  of  the  flow  over  the  type  of 
roughness  is  in  order. 

(7) 

The  data  of  Liu,  ct  ill.  were  obtained  on  transverse  square 
bar  roughness  with  a wavelength  of  4 in  low  speed  water  flow,  and 
their  results  indicate  an  equivalent  roughness  of  = 2 - 2.  5k,  fairly 
close  to  the  value  of  k^  = 2.  7 for  the  work  of  Betterman.  The  results 
of  Betterman  indicate  tliat  for  \ in  the  vicinity  of  4,  the  effectiveness 
is  extremely  sensitive  to  small  changes.  The  visualization  work  of 
Liu,  et  al.  revealed  that  in  the  vicinity  of  X = 4,  the  character  of  the 
flow  in  the  cavity  between  the  roughness  elements  changes  radically. 
For  X considerably  less  than  4,  the  flow  completely  bridges  the  cavity 

/ g \ 

(the  ”d"  type  flow  of  Perry,  et  al.  ' '),  creating  an  "open''  cavity, 
while  for  X much  larger  than  4,  the  flow  reattaches  to  the  cavity  floor, 
creating  two  distinct  separation  regions  within  the  cavity  (a  ''closed'’ 
cavityh  The  transition  from  one  type  of  cavity  to  the  other  occurs 
for  X somewhere  in  the  range  of  2 to  8.  This  explains  the  sudden 
shift  in  slope  of  the  correlation  results  of  Dirling  as  presented  in 
Figure  20  from  Mills  and  Courtney ^ ^ As  X increases  from  2,  the 
vortex  structure  within  the  cavity  is  stretched  and  pressure  communi- 
cation between  the  upstream  and  downstream  faces  is  decreased, 
resulting  in  increased  pressure  drag.  The  maximum  drag  is  reached 
very  near  the  value  of  \ at  which  rcattachment  occurs,  for  then  the 
pressure  drag  on  the  faces  has  reached  the  maximum  and  any  increase 
in  spacing  simply  produces  loss  force  per  unit  length,  and,  therefore. 
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less  drag  per  unit  length.  The  slope  change  occurs  at  a value  of  X 
slightly  above  4,  corresponding  to  the  change  between  open  and  closed 
cavity  flow, 

Charwat,  et  al.  have  investigated  higher  speed  flow  over 
two  dimensional  transverse  notches  and  conclude  that  for  supersonic 

i 

flow  over  large  cavities,  closed  cavity  flow  exists  for  X > 1^  and  open  ] 

cavity  flow  exists  for  X < 10.  The  observation  of  a crossover  point  ■ 

near  X = 4 for  subsonic  flow  is  also  made,  supporting  the  toncTusions  I 

of  Liu,  et  al.  j 

In  conclusion,  the  data  acquired  in  this  work  for  equ’librium  | 

rough  wall  flow  show  good  agreement  with  previously  published  data, 
when  correlated  with  respect  to  Reynolds  number  based  on  the  equivalent 
sand  grain  roughness  height.  The  equivalent  roughness,  however,  is 
found  to  be  significantly  altered  from  that  for  identical  roughness  in 
low  speed  flow.  These  data  support  the  claim  advanced  by  Mills  and 
Courtney , ^ ^ and  others  that  once  the  equivalent  sand  grain  rough- 
ness is  determined,  tlie  effect  of  a given  roughness  on  the  skin  friction 
and  velocity  may  readily  be  determined  (i.  e.  Figs.  18  and  19).  The 
determination  of  equivalent  sand  grain  roughness,  however,  is  not  a 
trivial  matter,  as  demonstrated  by  the  results  of  this  work.  In  par- 
ticular, the  extrapolation  of  any  low  speed  correlation  for  equivalent 
sand  grain  roughness,  such  as  that  of  Dirling,  to  higher  speed 

flows  must  be  viewed  as  highly  susceptible  to  very  large  errors. 

Fluctuation  Data 

Preliminary  analysis  of  the  data  during  the  acquisition  phase 
of  this  work  indicated  that  the  data  were  of  reasonable  magnitude  and 
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and  possessed  the  anticipated  profile  shape.  Later  complete  reduction 
and  extensive  analysis  revealed  that  some  of  the  data  were  obviously  : 

in  error,  but  it  was  not  possible  to  repeat  the  measurements.  As  a 
result,  fluctuation  data  for  some  of  the  wall  configurations  at  some 
axial  locations  were  discarded  and  are  not  included  in  this  discussion. 

The  fluctuation  magnitudes,  as  discussed  herein,  contain  rather 
large  error  bands- -probably  r 10-20^  on  a point  by  point  basis.  How- 
ever, when  considered  in  conjunction  with  data  at  other  points  through- 
out the  boundary  layer,  very  definite  trends  which  do  not  depend  on  the 
magnitude  accuracy  are  observed.  In  addition,  the  consistency  of  the 
observed  magnitudes  lends  considerable  credit  to  their  validity. 

Fluctuation  data  arc  presented  as  mass  flux  fluctuations 

>2  ~ ''  I — 

./m'^/m^  i,  total  temperature  fluctuations  ,/T  '*'/T  ;,  and  the 

' ^ 

associated  correlation  coefficient  R = m ' T ' / i I 'l 

’ mT^  ^t  ' W V t / 

for  all  wall  configurations  in  Figures  21,  22,  and  23.  No 

smooth  wall  data  were  obtained  at  x = 24,  4,  and  where  data  for 

a particular  configuration  have  been  omitted,  the  data  were  found 

to  be  in  error. 

Here,  as  for  the  mean  flow  data,  the  profiles  at  various 
axial  locations  were  compared  to  determine  where  the  profiles 
ceased  to  show  discernible  change.  The  boundary  layer  down- 
stream of  this  location  was  then  defined  as  an  equilibrium  or 
self  preserving  layer.  In  the  case  of  the  mean  flow  data,  the 
number  of  axial  locations  and  the  high  resolution  of  the  data  yielded 
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fairly  definitive  results.  For  the  fluctuation  data,  however,  data 
were  obtained  at  far  fewer  axial  stations  and  the  data  scatter 
was  much  greater.  As  a result,  the  axial  location  at  which 
the  profiles  become  self  similar  is  far  more  difficult  to  locate. 

Data  considerably  downstream  of  the  test  section  would  be 
required  to  verify  that  equilibrium  profiles  actually  are  estab- 
lished. 

Smooth  Wall 

The  smooth  wall  configuration  serves  as  a reference  for  deter- 
mining the  roughness  induced  effects,  and  as  such,  will  be  discussed 
first.  Examination  of  the  mass  flux  (m)  and  total  temperature  (T^) 
fluctuation  profiles  reveals  that  both  profiles  become  ''fuller'^  as  the 
axial  distance  increases  from  Z7.  4 to  39.  4 inches  from  the  throat  and 
then  remain  essentially  constant  downstream  to  x = 45.  4.  Thus,  al- 
though the  mean  flow  profiles,  as  discussed  above,  indicate  that  mean 
flow  equilibrium  on  the  smooth  wall  is  achieved  near  the  upstream  end 
of  the  test  rhombus,  the  fluctuation  properties  do  not  achieve  a self 
similar  or  equilibrium  state  until  some  distance  downstream,  between 
X = 33  and  x = 39,  some  14  to  20  inches  (18  to  256)  downstream  of  the 
test  section  boundary.  It  appears  that  the  fluctuation  levels  characteris- 
tic of  the  smooth  wall  are  suppressed  in  the  favorable  pressure  gradient 
flow  created  by  the  nozzle  and  gradually  relax  to  an  equilibrium  profile 
once  free  of  the  large  gradient.  The  data  of  Rose^^^^indicate  thatthe  fluc- 


tuation magnitude  in  a compressible  flow  is  significantly  affected  by  pres- 
sure gradients,  with  a favorable  gradient  causing  a decrease  in  magnitude. 
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Comparison  of  the  smootli  wall  da<-a  with  that  obtained  by 
Kistler/"^^^  L.adernian  and  Dernetriades^^^  ^ and  Owen  et  .il.  is 
made  in  Figure  24  (replottcd  from  Owen,  et  al.  ),  The  present  data 
are  observed  to  agree  quite  well  with  the  M = 1.  72  results  of  Kistler, 
in  opposition  to  liis  findings  that  all  (he  quantities  increase  with  iM.ich 
number  in  the  h'ach  number  range  of  1.7  to  4.  7.  Although  Owen, 
et  al.  attribute  the  differences  between  Kistler's  results  and  their  ciwn 
to  heat  transfer  effects,  the  present  data  do  not  substantiate  such  a 
claim,  for  in  spite  of  the  present  adiabatic  nozzle  wall,  the  data  do 
not  agree  with  the  trends  estabhsfied  by  Kistltir's  results. 

Smooth- to- Rough  T ransition 

Examination  of  the  smooth-to-rough  wall  configuration  data 
(Figs.  21,  22,  and  23)  reveals  that  self-similar  shapes  (within  the 
accuracy  of  the  d<ita)  are  attained  for  the  m fluctuation  profiles  up- 
stream of  X = 39  inches  for  all  roughness  sizes.  Equilibrium  profiles 
for  the  fliictuations  .ire  attained  somewhat  later,  between  x = 39.  4 
and  X = 45.4,  as  verified  by  the  available  x = 51. 4 data.  Following 
the  smootli-to-rough  wall  step  cliange,  then,  a definite  jiattern  is 
noted- -tHe  mean  flow  data  .ire  the  first  to  attain  new  equilibrium  pro- 
files, followed  by  the  mass  flux  fluctuation  data,  followed  in  turn  by 
the  total  temperature  fluctuation  data.  This  progression  is  a conse- 
quence of  the  manner  in  which  the  roughness  interacts  with  the-  boundary 
layer  to  produce  the  mass  flux  and  total  temper.ature  fluctuations.  The 
roughness  actually  induces  velocity  fluctuations  which  interact  with 
the  mean  flow  to  result  in  mass  flux  fluctuations  which  in  turn  interact 
with  the  mean  flow  and  velocity  fluctuations  to  creati;  total  temperature 
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fluctuations.  Thus,  before  the  fluctuations  caji  reach  equilibriiun  levels, 
the  mean  flow  must  be  in  equilibrium,  and  since  a finite  amount  of 
time  is  required  for  rh  fluctuations  to  respond  to  a change  in  me<m 
flow  conditions,  there  will  be  a time  (and  distance)  lag  prior  to  the 
establishment  of  in  equilibrium  profiles.  In  the  same  manner,  the 
fluctuations  are  dependent  upon  the  interaction  of  the  mean  flow 
and  rh  fluctuations,  and  so  will  reach  equilibrium  t'ven  later  in  time 
and  further  downstream. 

Roughness  effects  (the  increase  in  fluctuation  levels  with  re- 
spect to  those  over  a smooth  wall)  are  readily  apparent  at  x = 27.  4, 
remain  about  constant  in  the  inner  half  of  the  jirofile  but  increase  in 
the  outer  section  at  x = 33.  4,  and  then  decrease  somewhat  at  x = 45.  4. 
The  actual  rough  wall  fluctuation  magnitudes,  on  the  other  hand, 
develop  to  a self  similar  shape  prior  to  x = 39.  4.  The  apparent 
decrease  in  tlie  effdet  of  roughness  (smaller  amount  of  increase  due 
to  the  roughness),  then,  is  due  not  to  a decrease  in  the  rough  wall 
fluctuation  levels,  but  to  an  increase  in  the  smooth  wall  levels.  As 
discussed  previously,  the  presence  of  the  rough  wall  significantly 
alters  the  upstream  history  effects  of  the  mean  flow,  and  this,  in 
conjunction  with  its  incr(?ased  turbulence  production,  leads  to  an 
earlier  attainment  of  equilibrium.  Of  particular  interest  in  the  stream- 
wise  development  downstream  of  the  stef)  is  the  fact  that  the  m fluc- 


tuation level  at  the  wall  appears  to  first  decrease  upon  transition  to 
a rough  wall,  and  then  increase  along  with  the  rest  of  the  boundary 
layer.  The  total  temperature  data  indicate  no  such  decre.ise. 
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Careful  examination  of  the  data  reveals  further  insights  into 
the  manner  in  which  the  various  roughness  heights  affect  the  txirbu- 
lent  flow  structure.  Turning  first  to  the  rii  fluctuations,  consider 
Figure  21  in  which  these  data  are  presented  for  all  configurations. 

The  general  features  of  the  profiles  downstream  of  x = 27.  4 are  all 
similar,  but  those  at  x = 39.  4 illustrate  the  trends  as  well  as  any 
and  will  be  used  for  purposes  of  ilkistration.  The  largest  effect  is 
seen  to  be  caused  by  the  k = 0.  0125  inch  rough  wall  configuration 
which  also  exhibits  a conspicuous  maximum  some  y/6  = 0.  8 away 

from  the  wall  zis  opposed  to  the  more  ambiguous  maximum  exhibited 
by  the  other  configurations  located  at  distances  of  y/6  = 0.  5 for  the 

smooth  wall  and  k = 0.  050  inch  rough  wall  and  y/6  = 0.  8 for  the 

k = 0.  025  inch  rough  wall.  The  shape  of  the  k = 0.  012  5 inch  rough 
wall  mprofile  is  significantly  different  from  that  of  tlie  other  config- 
urations throxighout  the  inner  half  of  the  boundary  layer.  It  not  only 
has  a more  definite  maximum,  in  gener^^l,  as  noted  above,  but  it 
also  decreases  in  magnitude  very  r.ipidly  as  the  wall  is  approached- - 
so  rapidly,  in  fact,  that  near  the  wall  it  actually  falls  below  the  smooth 
wall  profile.  The  magnitude  also  decreases  rapidly  as  the  freestream 
flow  is  approached. 

The  k = 0.  025  inch  rough  wall  profile  data  also  exhibit  a peak 
magnitude  well  out  in  the  boundary  layer,  as  indicated  above.  This 
maximum,  however,  is  usually  less  pronounced  than  that  for  the 
k = 0.0125  inch  data,  with  the  magnitude  decreasing  less  rapidly  as 
the  wall  is  approached,  achieving  a value  somewhat  higher  than  the 
smooth  wall  configuration  profile  adjacent  to  the  wall.  Consideration 


-40- 


of  the  profile  between  the  maximum  magnitude  and  the  frcestream 
reveals  tliat  the  magnitude  drops  off  less  rapidly  than  the  k = 0.  01Z5 
inch  profile  as  the  freestream  is  approached  and  is  larger  than  that 
of  the  small  roughness  from  y/6  = 1.  2 to  the  freestream  at  y/6  =2.  0. 

Consider  now  tVie  k = 0.  050  incli  rough  wall  data;  it  is  observed 
that  its  maximum  occurs  much  nean^r  the  wall,  in  the  vicinity  of 
y/6  =0.  5,  the  magnitude  being  very  comparable  to  that  of  the 

k = 0.  02  5 inch  data.  In  fact,  the  behavior  of  the  profile  between  this 
maximum  magnitude  point  and  the  wall  is  virtually  identical  to  that 
of  the  k = 0.  025  inch  jjrofile.  doing  toward  the  freestream  flow,  the 
magnitude  first  drops  below  that  of  the  smaller  rougliness  configura- 
tions, and  then  assumes  virtually  the  same  value  as  the  k = 0.  0125 

>Jc 

inch  data  from  y/6  = 1,2  to  the  freestream. 

Thus,  although  the  data  for  actual  fluctuation  magnitudes  inay 

be  subject  to  considerable  error,  the  effect  of  roughness  sixe  is 

f \ 

definitely  disc€;rnible.  The  effect  of  a small  roughness  ( — — 7j  is 

w 

an  increase  in  the  m fluctuation  levels  throughout  the  bovindary  layer 
except  in  the  immediate  vicinity  of  the  wall,  wht;re  a slight  decrease 
may  occur.  A very  definite  maximnm  iTvagnitvide  j>eak  is  established 
well  displaced  toward  the  freenstream  compared  to  the  smooth  wall 
peak,  with  the  roughness  effect  decreasing  significantly  as  the  free- 
stream is  iipproached,  but  not  disappearing.  An  increase  in  the 
roughness  height  causes  an  increase  in  the  fluctuation  level  relative 
to  the  magnitude  peak  throughout  the  boundary  layer,  wdth  no  change 
in  peak  location.  As  the  roughness  size  continues  to  increase,  the 

j 

location  of  the  magnit\ide  peak  approaches  that  of  the  smooth  wall 
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magnitude  peak  and  the  fluctuation  level  between  this  peak  and  the 
freestream  decreases  somewhat  but  remains  above  the  smooth  wall 
level.  In  all  cases  the  rough  wall  freestream  fluctuation  levels  are 
about  double  the  smooth  wall  freestream  levels.  The  increased  le-vels 
throughout  the  flow  are  due  to  the  high  level  of  viTocity  fluctuation 
production  on  the  rough  walls. 

Turning  now  to  the  total  temperature  fluctuation  data  presented 
in  Figure  22,  it  is  apparent  that  significant  roughness  effects  an.- 
present  at  x = 27.  4 and  increase  in  the  outer  portion  of  the  boundary 
layer  prior  to  x = 33.  4,  then  decrease  somewhat  throughout  the 
boundary  layer  between  x = 33.  4 and  x = 39.  4.  Subsequent  to  x = 39.  4, 
little  change  occurs.  No  definite  roughness  size  effects  are  discern- 
ible in  the  data.  In  contrast  to  the  m fluctuation  data,  the  peak  magni- 
tude occurs  at  the  same  location  for  all  the  rough  walls  as  for  the 
smooth  wall. 

The  mass  flux-total  temperature  correlation  function  is  pre- 
sented in  Figure  23.  A rather  consistent  profile  is  seen  to  exist  at 
all  stations  for  all  the  svirface  conditions.  The  correlation  function 
starts  off  with  a small  negative  value  near  the  wall,  rapidly  becomes 
positive,  then  decreases  to  a rninimum  near  y/6  - 1.0  before  in- 

creasing again  and  approaching  a freestream  v.alue  of  about  0.  5-0.  7. 
The  smooth  wall  deviates  from  this  pattern  near  the  freestream 
to  assume  a negative  value  near  -0.  5 which  is  close  to  the  value  of 

/ -7  A \ 

-0.  7 foiuid  by  Owen,  etal.'  ^ at  M = 7 for  a non-adiabatic  boundary 


layer.  The  profile  shape  is  also  similar,  although  these  correlation 
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I'unctiona  arc  sii> nificantly  sn'ialler  througlicmt  the  body  of  the  tioundary 
layer. 

The  k = 0.  050  inch  rough  wall  data  generally  are  very  siniilar 
to  the  smooth  wall  data^  excej)t  near  the  bound. iry  layer  edge,  wl'ile 
the  smaller  roughnesses  maintain  somewhat  different  profiles,  just 
as  for  the  m fluctuation  profiles. 

Rough-to-Smooth  Transition 

Only  the  rough-to- smooth  step  change  fluctu.ition  data  pre- 
sented in  Figures  21-23  remain  to  be  analyzed.  Roth  the  m and  T^ 
fluctuation  profiles  at  x = 24.  4 are  observed  to  be  generally  similar 
to  those  for  the  k = 0.050  inch  smooth-to- rough  configuration  at 
X = 39  (fully  developed  or  self  similar  rough  wall  flow).  The  profiles 
are,  however,  concentrated  much  closer  to  the  wall,  probably  the 
result  of  the  strong  favor.ible  pressure  gradient  in  the  nozzle.  As 
indic^lted  earlier,  such  a gradient  tends  to  suppress  fhictuation  mag- 
nitudes somewhat.  The  correl.ation  functions  for  the  two  configurations 
are  very  similar. 

By  X = 27,  some  36  clownstream  of  the  step  change  in  roughness, 
the  character  of  both  fluctuation  profiles  lu-.s  changed  considerably, 
decreasing  in  magnitude  at  the  wall  <ind  increasing  in  magnitude  toward 
the  freestrearn.  This  growing  in  fhe  freestream  direction  accom- 
panied  by  develojjment  of  magnitude  peaks  near  y/6  = 0,  5 is  a reflec- 


tion of  the  adjustmt;nt  to  a zero  pressure  gradient  freestream,  while 
the  change  in  w.all  level  is  a reaction  to  the  shift  from  the  rough  to  the 
smooth  w.all.  The  correlation  fvinction  is  also  seen  to  respond  to  the 
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absence  of  the  favorable  pressure  gradient  by  (.'xp.uiding  in  tin*  free- 
stream  direction. 

The  profiles  continue  to  adjust  to  (he  changes  in  l)oth  the 
freestrearn  and  wall  conditions  downsti'eam  of  x = 27.  4.  A gradual 
increase  in  level  on  the  freestrearn  side  of  the  profiles  is  observed, 
accompanied  by  a shift  of  the  maximum  amplitude  peak  in  the  same 
direction.  Both  profiles  appear  to  be  close  to  equilibrium  at  x = 39.  4, 
but  continue  to  change  in  the  vicinity  of  the  freestrearn  between 
X = 39.  4 and  x = 45.  4.  A self  similar  profile  for  the  rrr  fluctuations 
appears  to  be  established  at  x = 45.  4,  but  the  profile  shows  evidence 
of  continued  evolvement  on  the  freestrearn  side.  Throughout  thi s 
adjustment  to  new  surface  conditions,  the  correlation  coefficient 
remains  essentially  unchanged  downstream  of  x = 27.  4. 

The  large  magnitude  of  the  rn  fluctuations  (larger  than  the 
k = 0.050  inch  smooth-to- rough  configuration)  is  due  to  the  locati<Jii 
of  the  rough  plate  in  the  nozzle.  The  same  roughness  height  was 
\ised  for  both  configurations,  but  the  rough-to- smooth  arrangem.ent 
placed  the  upstream  end  of  the  plate  in  a much  thinner  boundary  layer 
(6  ss  0.  2 inch  versus  6 0.  4 inch  for  the  smooth-to- rough  step 

change)-  The  effective  roughness  height  was  thus  much  greater,  well 
within  the  fully  rough  regiine.as  opposed  to  the  downstream  location 
where  it  was  at  the  very  lower  limit  of  that  regime.  The  turbulence 
production  caused  by  the  roughness  would  therefore  be  considerably 
larger  for  the  rough-to- smooth  model.  In  addition,  the  presence  of 
this  large  roughness  in  the  strong  pressure  gradient  section  of  the 
nozzle  would  significantly  pertvirb  the  mean  flow  and  exert  a larger 
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effect  on  the  upstream  boundary  history  than  wcnild  a similar  j)late 
in  the  zero  pressure  gradient  test  section. 

The  rough-to- smooth  configuration  fluctuation  data  <tre  also 
presented  in  Figure  24.  It  is  obvious  tliat  the  liistory  effects  are  very 
important  and  no  meaningful  comparison  of  nozzle  wall  data  may  be 
attempted  in  the  absence  of  more  detailed  knowledge  of  these  effects. 

In  summary,  the  fluctuating  flow  properties  were  found  to 
respond  to  changes  in  the  wall  roughness  sonaewliat  more  slowly  than 
did  the  mean  flow  properties.  This  plienomenon  is  due  to  the  mecha- 
nism by  which  tlie  fluctuating  properties  are  actually  influenced  by 
the  surface  condition.  Definite  roughness  size  effects  are  present  in 
the  rn  fluctuation  data  for  flow  over  a rough  wall,  but  none  were 
observed  for  the  T^  fluctuation  data.  Comparison  of  the  rough-to- 
smooth  step  change  data  with  the  smooth  wall  data  revealed  upstream 
history  effects,  just  as  mean  flow  data  did  earlier. 


-45- 


CONCLUSIONS 

An  cxpe rinientiil  invcjstifjation  of  tho  response  of  tlie  Hypersonic 
turbulent  boundary  layer  to  a step  cliangt;  in  surface  rtnighness  has 
been  performed  in  Leg  II  of  the  GALCIT  Hypersonic  V/ind  Tunnel. 

1 he  subject  boundary  layer  <leveloped  on  the  Hat  nozzle  wall  adjacent 
to  Mach  6 freestream  flow.  Both  mean  and  fluctuating  flow  data  were 
acquired  for  several  smooth-to- rouch  surface  step  changes  and  a 
single  rough-to- smooth  surface  step  change.  Use  of  a transverse 
square  bar  type  of  roughness  facilitated  direct  comparisons  witli  low 
speed  data.  Analysis  of  the  d^lta  from  this  program  has  yielded  the 
following  conclusions: 

The  establishment  of  new  equilibrium  mean  and  fluctuating 
flow  profiles  downstream  of  a step  change  in  surface  roughness  is 
accomplished  in  nearly  the  same  distance  (in  terms  of  boundary  layer 
thicknesses)  as  in  the  incompressible  case.  The  step  change  smooth- 
to-rough  configuration  boixndary  layer  attains  new  mean  flow  self 
similar  profiles  some  106  or  206  downstream  of  the  step,  while  the 
fluctuation  profiles  reach  this  state  some  14-166  or  28-326“^'  down- 
stream of  the  step.  These  distances  are  independent  of  roughness 

k u 

height  for  roughnesses  in  the  transitionally  rough  regime  ( ^<70). 

w 

The  step  roughness  change  rough-to- smooth  configuration  produced 
a somewhat  slower  adjustment  of  the  boundary  layer,  some  146  or 
26  6 for  mean  flow  equilibrium  and  some  20-256  or  40-50  6^'  for 

fluctuation  equilibrium. 


-46- 


2.  The  hypersonic  smooth  and  rougli  wall  equilibrium  profile 
velocity  data,  subsequent  to  a modified  Van  Driest  transf(jrir\at  ion 

to  "equivalent  incompressible"  form,  are  well  corx  elated  by  the 
-»  incompressible  composite  law  of  the  wall.  The  use  of  a virtual  oidgin 
for  the  effective  surface  is  necessary  to  achieve-  good  results  for  the. 
flow  over  the  surface  roughnt:ss  configurations,  just  as  for  incom- 
pressible flows.  The  modified  form  of  tlu;  Van  Driest  transformation 

4 

used  incorporated  a least  squares  quadratic  curve  fit  to  the  measured 
temperature  data  in  place  of  the  standard  Crocco  temperature  profile. 

3.  The  above  correlation  is  also  found  to  be  applicable  in  tin: 
highly  non-equilibrium  flow  region  immediately  downstream  of  the 
step  change  in  roughness,  provided  the  rougliness  is  not  too  large. 

A roughness  height  of  0.050  inch  (k^  =.  45)  for  the  smooth-to- rough  step 
surface  change  was  found  to  be  too  large  in  tliis  sense,  but  the  same 
height  used  in  the  rough-to- smooth  step  change  configuration  was  not. 

4.  The  effects  of  surface  roughness  on  c?quilibriuin  or  self- 
preserving boundary  layers  were  found  to  exhibit  excellent  agreement 
with  previously  p\iblished  data,  when  correlated  versus  equivalent 
sand  grain  roughness.  Wlien  viewed  in  this  manner,  with  the  e(juiva- 
lent  sand  grain  roughne.ss  a unique  parameter  for  a given  type  vif 
roughness  in  a given  freestream  velocity,  the  effect  of  a change  in 
roughness  height  on  the  skin  friction  and  nondimensional  velocity 
(u/u.^)  profiles  is  found  to  be  a universal  function. 

5.  Kxisting  procedures  (based  on  incompressible  data)  for 
calculating  the  equiv.ilcmt  sand  grain  roughness  of  a given  roughness 


configuration  were  found  to  he  inapplicable  in  the  present  case,  due 
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to  velocity  induced  changes  in  the  flow  over  the  roughness.  In  light 
ol  this  finding,  the  extrapolation  of  any  low  speed  rougluiess  corre- 
lation to  higher  velocities  must  be  viewed  very  critically  with  the 
knowledge  that  £iny  results  of  such  a procedure  are  very  susceptible 
to  gross  errors.  A concentrated  effort  to  determine  tlie  effects  of 
compressibility  and  velocity  on  such  correlations  is  definitely  needed. 

6.  The  investigation  of  the  rough- to- smooth  step  change 
configuration  revealed  the  existence  of  very  significant  pressure  and 
temperature  history  effects  throughout  the  boundary  layer.  The 
existence  of  these  effects  affords  an  explanation  of  apparent  incon- 
sistency of  data  obtained  on  wind  tunnel  nozzle  walls  with  that  obtained 
on  flat  plates  supported  in  the  tunnel  freestream.  The  presence  of 
the  rough  surface  on  the  nozzle  wall  within  the  expansion  region  was 
sufficient  to  alter  the  history  effects,  but  was  not  able  to  eliminate 
them.  The  unknown  nature  and  extent  of  such  history  effects  casts 
serious  doubts  on  the  validity  of  simulating  boundary  layers  on  free- 
stream surfaces  by  use  of  nozzle  wall  boundary  layers.  Considerably 
more  research  is  needed  in  this  area  to  delineate  the  extent  of  these 
history  effects  and  investigate  the  possibility  of  destroying  them. 
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